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i — i Abstract 

Q_i Equations of magneto-gasdynamics in the natural curvilinear system of co- 

ordinates where trajectories and magnetic lines play a role of coordinate 
curves are reduced to the nonlinear vector wave equation coupled with the 
incompressibility condition in the form of the generalized Cauchy integral. 
The symmetry group of obtained equation, equivalence transformation, and 
t-h group classification with respect to the constitutive equation are calculated. 

New exact solutions with functional arbitrariness describing non- stationary 
incompressible flows with constant total pressure are given by explicit formu- 
lae. The corresponding magnetic surfaces have the shape of deformed nested 
cylinders, tori, or knotted tubes. 
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^ Introduction 

The topological structure of the magnetic field plays a significant role 
in applied problems (magnetic traps for plasma confinement, MHD-dynamo, 
etc.) and related mathematical theories [H El E]- Most of these theories re- 
quire the magnetic field to be determined in a finite area of the 3D space, 
and satisfy natural conditions over the boundary of the area. However, the 
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set of known examples of exact solutions of magneto-hydrodynamics equa- 
tions (MHD) that describe such topologically nontrivial configurations is very 
poor. The main problem is that the existing methods (e.g. symmetry anal- 
ysis) do not allow finding exact solutions with prescribed geometrical prop- 
erties of the governed flow. The goal of this paper is to develop a suitable 
framework for the description of MHD flows with non-trivial topology, and 
to construct new classes of exact solutions that possess functional arbitrari- 
ness and govern non-stationary plasma flows with knotted configurations of 
magnetic lines and magnetic surfaces. 

The main idea of the paper is to "hide" the geometry of the flow into 
the special choice of a curvilinear system of coordinates. Namely, the system 
of coordinates is chosen such that trajectories and magnetic lines of the 
flow form two families of coordinate curves. This approach allows one to 
separate the topological structure of the magnetic field and the evolution 
of the flow. The obtained system of coordinates is referred to as "natural" 
by the analogy with the natural coordinates in classical mechanics. In the 
natural coordinates the system of MHD equations is reduced to a non-linear 
vector wave equation and to the incompressibility condition in the form of 
generalized Cauchy integral. The symmetry group of the obtained system 
of equations, equivalence transformations, and a group classification with 
respect to the constitutive equation of fluid are calculated. 

The obtained system of equations allows construction of a class of exact 
solutions describing non-stationary incompressible MHD flows with constant 
total pressure. The complete description of this class of solutions requires 
separation of variables in a certain scalar equation with the subsequent so- 
lution of an overdetermined system of partially differential equations. This 
procedure in particular cases gives new exact solutions of MHD equations 
possessing a functional arbitrariness up to one function of three variables, 
two functions of two variables and one function of one variable. The so- 
lutions describe non-stationary plasma flows with cylindrical, toroidal, or 
knotted magnetic surfaces. 

The stationary case was investigated in the similar manner in Re- 
markably, that in the analogous natural system of coordinates where stream- 
lines and magnetic lines form two families of coordinate curves, MHD equa- 
tions were also reduced to a vector wave equation subjected to a geometrical 
constraint. The class of solutions with constant total pressure was also de- 
scribed. It was shown, that "Maxwellian" surfaces that bear magnetic lines 
and streamlines of the flow belong to the class of translational surfaces ob- 
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tained by sliding one 3D curve along another. Similar system of coordinates 
for stationary MHD equations was used by C. Rogers and W. K. Schief [5] 
where the relation between solutions in which the total pressure is constant 
on Maxwellian surfaces with the integrable Pohlmeyer-Lund-Regge model 
was discovered. The same method was applied to the solution of Gilbarg 
problem [6] analogue: to which extent the geometry of the flow (streamlines 
and magnetic lines pattern) determines the flow, and for the construction of 
one-parameter family of exact solutions describing an equilibria state with 
toric geometry [Tj. In the non-stationary case the similar system of curvi- 
linear coordinates constructed on trajectories and magnetic lines of the flow 
was used by J. D. Gibbon and D. D. Holm [5] for the analysis of alignment 
dynamics of non-stationary MHD flow using a quaternionic approach. 

1. Preliminary calculations 

Ideal magneto-gasdynamics equations have the following dimensionless 
form: 

p t + u ■ Vp + p div u = 0, (la) 

p(u t + (u • V)u) + B x rotB + Vp = 0, (lb) 

B ( = rot(uxB), (lc) 
divB = 0. (Id) 

Here u is the velocity vector, B is the magnetic field, p is the density, and p 
is the pressure. System of equations ([I]) in the case of compressible fluid is 
extended by the state equation p = F(p, S) with entropy S that is conserved 
along particles' trajectories: 

S t + (u • V)S = 0. 

For the incompressible fluid p = const the state equation is not required. 
Pressure p acts as an unknown function. 

The system of MHD equations ([I]) is transformed as follows. By using 
the vector identity 

B x rotB = ^V|B| 2 - (B ■ V)B, 
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and introducing the total pressure P = p + ||B| 2 the momentum equation 
(lb) is brought to the form 

p(u t + (u ■ V)u) - (B ■ V)B + VP = 0. 

With the use of an identity 

rot (u x B) = u div B — B div u + (B ■ V)u - (u • V)B, 



and equations (la), (Id) equation (lc) is modified as follows: 



B 

P 



(pt + u- Vp) + (B ■ V)u - (u ■ V)B. 



The latter can be conveniently written as 



(2) 



where notations 



a + u ■ v, d s 



B 



for differentiations along trajectories and magnetic lines respectively are in- 
troduced. Note, that equation ^ is equivalent to the commutativity of 
vector fields D t and D s in IR 4 (t,x). 

Finally, with the use of the modified magnetic field b = p _1 B equations 
([!]) are written as follows: 



p t + u • Vp + p div u = 0, 

p(D t u-D 3 (pb)) +VP = 0, 

D t b = D s u, 
div (pb) = 0. 



(3a) 

(3b) 

(3c) 
(3d) 



This form of MHD equations will be used for construction of the curvilinear 
coordinate system in the following sections. 
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2. Natural curvilinear system of coordinates 

A curvilinear coordinate system (£,£), determined for each fixed t by 
the diffeomorphism x = j(t, £) is introduced over an open subset of the 
Euclidean space IR 4 (t,x). Basic vectors of the coordinate system are defined 
at each point as 

e o = (1, 7t)> e i = (°> 7e)> * = 2 ' 3 - 

As before, the lower indices denote the corresponding partial derivatives. 
Commutation relation ^ for vector fields D t and D s allows choosing these 
fields as basic ones [9]: 

e = (l,u), ei = (0,b). 

The two remaining vector fields e 2 and e 3 are not fixed. Thus, coordinate 
lines of t-family coincide with particles' trajectories; Coordinate lines of £ l - 
family trace the magnetic lines. At that, relation ^ is equivalent to the 
equality of mixed derivatives of function 7: 

d dj <9 2 7 d dj 
d£ l dt dtd^ 1 dt d^ 1 

In what follows the notation £° = t, and the Einstein convention of sum- 
mation with respect to the repeating indices are adopted. At that, if the 
summation is performed with respect to a Greek index, it is supposed to 
take values from to 3; The repeated Roman index runs from 1 to 3. Con- 
travariant components of vectors u = u a e a and b = b a e a in the curvilinear 
coordinates take the following simple form: 

u° = 1, u l = 0, i = 1, 2, 3; b° = 0, b 1 = 1, b 2 = b 3 = 0. (4) 

The constructed curvilinear system of coordinates will be referred to as 
'natural' by the analogy with natural coordinates in classical mechanics. The 
natural system of coordinates can also be regarded as the special choice of 
Lagrangian coordinates. The arbitrariness in the choice of Lagrangian co- 
ordinates £ is used for "straightening" of magnetic field lines. Thus, the 
complexity of the magnetic field is hidden in the proper choice of the coordi- 
nate system. The initial magnetic field specifies the coordinate system in the 
initial state. Time evolution of the magnetic field automatically follows from 
the evolution of the coordinate system due to the frozenness of the magnetic 
field into the flow. 
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3. Equations in natural coordinates 

In this section equations ^ will be written in the natural coordinates with 
the use of formulae given in Appendix Appendix A By virtue of expression 
(El) one has 



div u 



U A 



u p T\ 



' 0i y/g dt' 

Here f i denotes the covariant derivative of function f with respect to i- 



th variable; g is the determinant of the metric tensor; T^g are Christoffel 

(5) 



symbols of second kind. Equation (|3a|) is transformed as follows 

ln(pv^) =0. 



dp P d^/g 
dt y/g dt 







dt 



In the same way, equation (pdl can be written as 



(P&*),< 



dp b 



- + pb a T 



at 



dp 
dp 



P d^g 



+ 



d 



Equations (pj) and ^ lead to the generalized Cauchy integral 

P^g = f(i\e) 

with an arbitrary function /. 



T m(/Vs)) =0- (6) 



(7) 



It is remained to rewrite the momentum equation (3b). The analogous 
computation of covariant derivatives taking into account formulae (|4]) gives 



Note the following identities: 
loo*- Qt , 



0. 



1,2,3. 



(8) 



r il ej 



dei 



^g 



ik 



Multiplication of equations (|8j) on basic vectors and summation with re- 
spect to i in view of the preceding identities gives 



P 



9 2 7 d 



dj 



dj 



VcP = 0. 



Hereafter symbol "T" denotes the transpose of the matrix; (d£/d"y) is Jaco- 
bian matrix. Finally, equations for the sought mapping x = *y(t, £) take the 
following form: 

pVg = f(e,e), v^=det^. 

In the case of the incompressible fluid (p = p ) function P is treated as 
an unknown function. For the compressible fluid with the state equation 
p = F(p, S) relations ^ are extended by the expression for the complete 
pressure in terms of density and entropy as 

(10) 

with an arbitrary function S(£). Thus, in the natural system of coordinates 
original MHD equations ([I]) are reduced to the nonlinear vector wave equation 
and to the incompressibility condition in the form of the generalized Cauchy 
integral. 



4. The physical picture of the flow and initial conditions 



Each solution x = 7(i,£) of the system (J9J), (10) provides an explicit 
description of particles' trajectories and magnetic lines of the governed fluid 
flow. 

Particle trajectory originating at time t — at the position x = 7(0, £ ) 
is given explicitly by the parametric formula x(t) = "f(t,$, ). Magnetic lines 
are parameterized by coordinate £ 1 . At the arbitrary moment of time to the 
parametrization of a magnetic line has the form x(s) = 7(to> s, £q, £o)- At 
that, coordinates £q and £q "distinguish" different magnetic lines. 

Note, that the initial function 7(0, £) = 7q(£) can not be taken arbitrarily. 
Indeed, let the initial density distribution p(0,x) = po(x) be given. Then, 
the dependence 7 (£) should satisfy the equation 

fl)(7o(0)det(^) =m 2 ,e) (11) 
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which is equivalent to the incompressibility of the initial magnetic field stated 



by equation (Id) 



Suppose the divergent-free magnetic field B (x), density po( x ), and initial 
velocity field u (x) of fluid be given at initial time t = 0. Hence, the magnetic 
field uniquely covers the area of the flow such that only one magnetic line 
passes through each point. Construction of the initial data for system ^ is 
performed as follows. 

With the vector field b (x) = p 1 (x)B (x) one solves the system of ordi- 
nary differential equations 

dx u f \ 
dfi = b oW- 

The general solution of this system contains three arbitrary constants. One of 
these constants corresponds to the trivial admissible translational symmetry 
^ — y ^ + a and is not important. The two remaining constants are denoted 
as £ 2 and £ 3 . As a result, this gives a diffeomorphism x = 7 (£) (at least, 
for smooth vector field b ). Moreover, due to the divergence- freeness of 
the initial vector field B (x) this diffeomorphism automatically satisfies the 



incompressibility condition (11). 



For the given initial velocity field Uo(x), the vector field 

7i«)=uo(7o«)) 

is constructed. This finally gives the initial conditions for the system of 
equations ^ in the form 

7((U)=7o(£), ^(0,£)=7 a (£)- 

In case of incompressible fluid these data are extended by specification of the 
initial pressure: 

m£)=Po(7o(£)) + ^|£o(7o(£)) 



5. Symmetries and equivalence transformations 

System of equations ^ is related to the initial MHD system in Cartesian 
coordinates by the non-local transformation and partial integration. Hence, 
the admissible group of symmetries is not necessarily conserves. In this 
section symmetry properties of equation ([9j are investigated. 
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Equations (|9| admit the following equivalence transformation 



d(P,P) T 



(12) 



Here d(-, -)/d(-, ■) denotes the corresponding Jacobian. Since / 7^ 0, by the 

2 3 

suitable choice of functions £ (£ 2 ,£ 3 ) and £ (P,P) one can make f — 1. In 
addition to the equivalence transformation (12), system Em always admits 
the "reparametrization" symmetry 



1. (13) 



—2 —3 

with arbitrary smooth functions a, £ and £ . 



Symmetries of system (|9j) are calculated separately for cases of compress- 
ible and incompressible fluids. In the incompressible case the density is nor- 
malized to unity: p — 1. In order to have an explicit expression of the time 
derivative it is convenient to rewrite system ^ in the following form: 



dt 2 



d 2 j dP 
+ 



dj dj 

X 



dp 2 dp \dp «9£ 
dP ( dj 



+ 



OP 

dp 



<?7 dj 

dp x W 



+ 



dp \dp dp 
dj dj dy 
dp' dp' d~P 



(14) 



1. 



Here [-,-,•] denotes the mixed product of three vectors; x is the vector prod- 
uct in IR 3 . 

The direct computation of the symmetry group of system ( 14 ) using the 
standard algorithm [TQl E] gives the following result. 



Theorem 1. The symmetry group of system (14) is a semidirect product 
of the reparametrization group (13) and the generalized Galilean group. The 
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latter is generated by the following infinitesimal operators: 



T 



dt, 



7*- 1 



1,2,3, 



S 2 = 3£ 2 <% + 3£ 3 <% + 2f'<9y + 4Pd P , 
Gi{a\t)) = d L {t)d^ -7 i a i (t)9 P , i 
K = p(t)d P . 



1,2,3, 



Here the summation is taken only with respect to j,k = 1,2,3; a % {t) and 
(3(t) are arbitrary smooth functions; is the Levi-Civita symbol; upper dot 
denotes the t-derivative. 

Note that the admissible group, except for its "reparametrization" part 
(13), is inherited from the original MHD equations ([I]) in Eulerian coordi- 
nates. Operator T generates a time shift, operators Ri give the group of 
rotations in M 3 (7), S\ and S 2 generate dilatations, operator K claims that 
pressure can be modified by addition of an arbitrary function of time. The 
most non-trivial transformation is the generalized Galilean translation gen- 
erated by operators Gi(a l {t)) 



^ = 7 + CK (t), P = P - ~i ■ ix{t) - -cc{t) ■ cc{t) 

with an arbitrary time- dependent vector-function a.(t). 

In the case of compressible fluid with the constitutive equation p 
F(p, S) relations ^ can be conveniently written as follows: 



<9 2 7 



<?7 



d_ 

dt 2 df 1 
dP 

+ w \de 

dj dy dy 

W'de'de 



+ 



X 



dP ( 


^7 


~W\ 


d~e 


de) ' 


= o, 


dp 
' di 


= h 



+ 



dP 



d*y 

d~e j 1 a? 



dj dj 

d~e x d~e 



(15) 



h(p,p) 



dp 

di' 



v + - 2 p 



dj 



de 



Here h(p, p) is a square of the speed of sound h(p, p) 
functions in system (15) are 7, p and p. 



dF I 



dp \S=S(p,p) 



. Unknown 
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Computation of the symmetry group of system (15) requires the group 
classification with respect to the constitutive equation represented by func- 
tion h(p,p). The system admits the following equivalence transformations: 

t = a% f = e, f = P 3 e, ? = P 3 e, 7=/? 2 7, p=<*-*p, 

p = a - 8 f3 4 p + K } h(p,p) = a- 2 (3 4 h(a s (3- 4 {p - K),a 6 p) 

(16) 

with arbitrary constants a, (3, and K. The result of calculation of the sym- 
metry group is the following. 



Theorem 2. The kernel of admissible group of equations (15) is a semidi- 
rect sum of the infinite-dimensional reparametrization group (13) and the 
extended Galilean group generated by the following infinitesimal operators: 

T = d t , 

T% — — d>yi ? % — X ? 2 ? 3 j 

Ri = Sijk^dyk, 2 = 1,2,3, 

Gi = td^i, i = 1,2, 3. 

As before, the summation is taken only with respect to j,k = 1,2,3; 
is the Levi-Civita symbol. This group is admitted for an arbitrary function 
h(p,p). For special functions h the kernel of admissible group is extended by 
the family of operators 

Y = atdt + (-Act + 2c 2 )£ 1 d e + c 2 £ 2 <% + c 2 £ 3 «% 

+ 2c a ya y i + (c 3 + (4c 2 - 8 Cl )p)d p + (4c 2 - 10 Cl )pd p 

provided the following equation is identically satisfied for p and p 

2cx(h + Aph p + 5ph p ) - Ac 2 (ph p + ph p ) - c 3 h p = 0. (17) 

Here Ci are arbitrary constants. 

Group classification with respect to h(p, p) is performed in the same way 
as in classical gas dynamics [TU] . The result is summarized in Table [l] Here 
/, F, and g are arbitrary functions, k is an arbitrary real parameter. The 
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Table 1: The result of group classification 





ti{p, p) 


r\Pi a ) 


Vjrl (J Lip cJvLcIlolUIl 


1 





F(S) 


Yi, Y 2 , Y 3 


2 


p k , k^-1 


^- lP k+1 + F(S) 

K-\- 1 ' v ' 


2kY 1 + (1 + 5k)Y 2 , Y 3 


3 


p- 1 


\np + F(S) 


Y t + 2Y 2 ,Y 3 


4 


kp/ p, k 7^ 


F{S)p k 


Y X ,Y 2 


5 


P k f(pp k ~ 1 ) 


g(pp k - l )F(S) 


2kY 1 + (1 + Ak)Y 2 


6 


f(pe p )/ P 


g{peP)+F{S) 


Y 1 + 2Y 2 + 2Y 3 


7 


f(p) 


9(p)+F(S) 


Y 3 


8 


f(p)/p 


g{pHS)) 


Y 1 + 2Y 2 


9 


h(p, p) 


F(p,S) 





first column represents the order number of the group extension. The sec- 
ond column gives the particular form of function h(p, p) normalized by the 
equivalence transformations ( Jl6) ). In the third column the corresponding 
constitutive equation p = F(p, S) is specified. The last column of the Table 
provides infinitesimal operators of transformations, that extend the kernel of 
the admissible group for the particular constitutive equation. The following 
notations for the extension operators are used: 

Yi = td t - 4^% + 2yS^ - 8pd p - 10pd p , 

Y 2 = 2^d e + + (% 3 + Apd p + 4pd p , Y 3 = d p . 

Operators Y\ and Y 2 correspond to dilatations, operator Y 3 generates the 
shift of pressure. Details of the group classification can be found in Appendix 
Appendix B 



6. Flows with constant total pressure 

This section is devoted to the construction of a special class of solutions 



to equations ( 14 ) of incompressible infinitely conducting fluid. Suppose that 



the total pressure is constant over the whole area occupied by the flow: 

1 |2 

P = p + — |B| = const. 
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From the first of equations ( 14 ) it follows that 



dt 2 



0. 



Integration of this equation yields the expression for function ~f(t,£): 
With this representation the incompressibility condition (the second of equa- 



tions (14)) must be satisfied: 



-<7i + Tj 



<T 2 + T 2 ) X (<T 3 + T 3 



1. 



(19) 



Hereafter the lower index l i' denotes the derivative with respect to z-th ar- 
gument of functions cr and r. In the expanded form scalar equation (19) 



contains 48 additives. Each of them has the form of a product of two multi- 
pliers depending on different variables: t — and t + C, 1 . Thus, it is required 



to separate variables in equation (19) and then to solve the obtained overde- 



termined systems of nonlinear partially differential equations for functions 
r and cr. Not pretending to solve this complicated problem in the general 
form, two nontrivial examples of exact solutions are constructed below. 



7. Stationary field-aligned flows 

Suppose that a = over the whole area of the flow. The solution is any 
vector field r with a unit Jacobian: 

7(t, o = r(t + e, e, a, det = 1. (20) 

In this solution velocity vector u and magnetic field B are collinear. The 
solution describes stationary flow of plasma aligned with the magnetic field 
known as Chandrasekhar's solutions [T2] . 

8. Non-stationary flows 

In this special case it is supposed that vector field cr depends only on one 
variable t — C, 1 : 

1 = (7 (t-e)+r(t+e,e,e). m 
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Equation ( 19 ) is reduced to the following one: 

- cr' ■ (r 2 x r 3 ) = 1 - n ■ (r 2 x r 3 ). 



(22) 



Prime denotes the derivative of vector field cr with respect to its only argu- 
ment. 

Henceforth notations 77 1 , r] 2 , and T7 3 for vectors of standard Cartesian 
basis in IR 3 are adopted. Vector field r has the following decomposition: 
r = 5^i=i rlr 7 l - The solution of equation (22) depends on the dimension 
dim{cr} of the linear space spanned by the vector cr'(t — £ x ) as its argument 
varies. 

dim{cr} = 1. Without loss of generality it is supposed that vector cr' is 



equal to T7 1 . Separation of variables in equation (22) gives 



tr , = rj L , (n - ry 1 ) • (t 2 x T 3 ) = 1. 



Integration of the first equation with respect to t — £ x accurate to the insuf- 
ficient constant additive vector yields 

a={t-ew. 

It is convenient to represent vector r as r = (t + ^^ + T. Then the solution 
is reduced to (tilde is dropped) 



Now we see that accurate to the Galilean translation along 77 1 this solution 



coincides with the stationary solution given by (20). 



dim{cr} = 2. Separation of variables in equation (22) gives 



U(t - eW + V 2 , V 1 ■ (r 2 x r 3 ) 



[Ti - rj ) ■ [t 2 x r 3 ; 



1. (23) 



Here U is a non-constant function. Integration of the first equation of (23) 
accurate to the insufficient additive constant vector leads to 
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Here u'(a) = U(a). The second of equations (23) implies that components 
t 2 and r 3 of vector r are functionally dependent as functions of £ 2 and £ 3 : 



r 2 = r 2 (t + e 1 ,A(t + e 1 ,e 2 ,a): 



It is convenient to change vector r as r — (t + C, 1 ) 1 ! 2 — > t. At that, vector 7 
acquires an additive 2£r/ 2 , which can be zeroed by the Galilean translation. 
Finally, vector 7 is obtained in the following form: 



7 = (u(t - e) +r\t + £\ e, £ 3 )y + r 2 (t + £\ A(f + £\ £ 2 , £ 3 ))*7 2 

+ r 3 (t + e 1 ,A(t + e 1 ,e 2 ,a)^ 3 ) 



7T rr 



7T 



A? Aa 



1. 



(24) 
(25) 



Scalar equation (25) states that the product of two Jacobians must be equal 



to 1. Note, that without normalization (12) the right-hand side of this equa- 
tion is an arbitrary function /(£ 2 ,£ 3 )- Solution (24) contains one arbitrary 



function of three arguments (either of functions r 1 or A), two arbitrary func- 
tions of two arguments (r 2 and r 3 ), and one arbitrary function of one argu- 
ment (function u). Equation (25) can be treated as a linear equation for r 1 
provided all remaining arbitrary functions are specified. 



dim{cr} = 3. Separation of variables in equation (22) gives 

a' = U\t - eW + U 2 (t - eW + l\ V 1 ■ (r 2 x r 3 ) = 0, 

rf ■ (r 2 x r 3 ) = 



77 ■ 

(Tl 



Vi) ■ (r 2 x t 3 ) = 1. (26) 



Here functions U , U , and 1 are supposed to be linearly independent as 



functions of t—^ 1 . Integration of the first of equation (26 ) gives the expression 
for vector cr: 

a = u \t - + « 2 (t - ew +(t- ew- 



The second and the third equations of (26) are compatible with the last 
equation in (26) only if r 3 does not depend on £ 2 and £ 3 : t 3 = r^it + t; 1 ). As 



before, accurate to the Galilean translation, one obtains the solution: 
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+ (u 2 ( t - e) +r 2 (t+ a) v 2 + At + ew, ^ 
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(28) 



This solution contains one arbitrary function of three arguments (either of 
r 1 and r 2 ), and three arbitrary functions of one argument (functions r 3 , u 1 
and u 2 ). Scalar equation (28) can be treated as a linear equation for r 1 
that also contains an arbitrariness in one function of two arguments. Both 
solutions can be extended by the application of the admissible symmetry 
transformations specified in Theorem [T] 

Solutions (24) and (27) can be treated as non-stationary disturbances of 



the field-aligned flow (20). The overall geometry of the flow is governed by 



the vector field r, while functions u and u % can be specified arbitrarily and 
independently providing perturbations of the stationary flow. 

9. Some particular examples of solutions 



In this section examples of plasma flows governed by solutions (24) and 



(27) for particular choices of arbitrary functions are constructed. 



The first solution under consideration is (27). For simplicity the normal- 

mation 

dr 1 dr 2 



ization (12) is not used, so that equation (28) reads 

dr 1 dr 2 



Q£2 Q£3 Q£3 fl£2 



(29) 



with an arbitrary function /. 

On the first step the basic field-aligned stationary solution is constructed 
by setting u l = 0, u 2 = 0. Functions r 1 and r 3 are taken as 



T 



a(n) + a(fjL) (A(e, e) cos <p(ji) + B(e, e) sin (p(fj)) . 

r 2 = bo*) - m {A(e,e) am <p(n) - B(e, e) «» , 



(30) 



Equation ( 29 ) after separation of variables yields 
dr 3 



d/x 



d(A,B) 



(31) 



The first equation in (31) determines functions r 3 in terms of functions a 
and /3; The second equation specifies functions /(^ 2 ,^ 3 ). 
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Solution (30) can be treated as follows. Functions A and B can be chosen 



arbitrarily to select a family of curves in the plane spanned by vectors T7 1 , 
1 For example, the choice 



77 



A = k^ 3 cos £ 2 , B = k 2 £ 3 sin £ 2 



specifies a family of nested ellipses parameterized by £ 2 with semi-axes k^ 3 
and k 2 £, 3 . This planar picture of curves propagates along rf direction being 
rotated on angle ip(fj), stretched and shifted in rf and rf directions as spec- 
ified by functions a(/i), P(fj), a(/i), and b(fi). Since a and are non-zero 
functions, dependence t 3 (/i) is monotonic, therefore variable \x (or C, 1 ) can be 
treated as a parameter along rf. Magnetic surfaces are obtained by fixing 
either £ 2 or £ 3 . In general, this procedure gives a set of nested deformed 
cylindrical surfaces for the magnetic surfaces in the basic stationary state. 

On the second step one chooses non-zero functions w 1 (t — u 2 (t — £ 1 ) 
to obtain the non-stationary disturbance of the constructed basic stationary 
state. As coordinate plays a role of a parameter along r] 3 , the disturbance 
propagates along rf as time grows. 



Summing up, solution (27) describes non-stationary plasma jet flows 



stretched along one spatial direction. 



Now solution (24) is observed. Again, the normalization (12) is not per- 



formed, so that equation (25) reads 



dr 2 dr 3 dr 2 dr" 



dX dr 1 dX dr 1 



(32) 



dC dX dX J \dZ 2 d£ 3 d£ 3 d£ 2 . 

On the first step the basic stationary solution is constructed by choosing 
u — 0. Arbitrary functions in p4j) are selected as follows: 



r 1 = (aO + B(e, £ 3 ) sin^/i) + A(£ 2 , £ 3 )) , 



A = yjb(ji) + 5(e 2 ,e) cos(^(/i) + A(e^ 3 )) 
t 2 = X cos kfi, t 3 = Asin/c/i, / u = t + £ 1 . 



(33) 



Substitution of (33) into (32) gives 

f(e,e) = 



kB d(A, B) 



which can be treated as the definition of function /. Solution (33) describes 



stationary flow with nested toroidal (or knotted) magnetic surfaces provided 



17 




Figure 1: The magnetic surface £ d = 1 with the magnetic line in the shape of the trefoil 



knot described by solutions ( 33 ) , ( 34 ) 



the inequality |£>(£ 2 ,£ 3 )| < b(fj) is satisfied. At that, variable \i grows along 
the central curve of the torus; ^p(fi) is responsible for the winding of the 
magnetic lines around the torus, functions A and B determine the shape of 
the section of the toroidal surface, functions a and b describe the shift of the 
section of the torus relative to its central curve. 



Suitable choice of arbitrary functions in solution (33) yields various pic- 
tures of magnetic surfaces. For example, the choice 



A = £ 2 , B = £ 3 , <p = 3/i, a = 0, b 



k 



(34) 



generates the toroidal magnetic surface in Figure [T£i. All magnetic lines on 
this surface are closed curves topologically equivalent to the trefoil knot. 
The choice 



A = e, s = e, ¥>=3//, 

6 = 3 + cos3/i, k = 2 



sin 3/x, 



(35) 



leads to the magnetic surfaces in the shape of the trefoil knot depicted in 
figure [T]d. One can easily generalize these examples to more topologically 
complicated pictures of plasma flows. 
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Figure 2: The magnetic surface £ 3 = 1 in the shape of the trefoil knot described by 
solutions (33), (35). 



On the second step one can add a non-stationary disturbance of con- 
structed stationary pictures by choosing a non-zero function u. The distur- 
bance propagate along the central curve of the magnetic tubes as time grows. 
Periodic functions u generate time-periodic disturbances. 

As a conclusion, solution (24) describes various plasma flows with torus- 
or knot-shaped magnetic surfaces. 
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Appendix A. Some useful formulae 

In this section some useful facts from differential geometry are reminded. 
Along with the basic vectors e Q defined in section [2j one can use vectors of 
the cobasis e a . To this end the smooth inverse mapping £ = 7~ 1 (t, x) is 
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utilized. The cobasic vectors are determined as 



e° = (l,0), e l = &,V x e), 2 = 1,2,3. 

Co- and contravariant components of the metric tensor are calculated ac- 
cording to 

ad a B 

g a p = e a -ep, g p = e • e p . 

The Christoffel symbols of the second kind are determined by the expansion 
of derivatives of basic vectors: 

de a _ p dea _ v b 
Note, that the special form of basic vectors e a yields 

The expressions for Christoffel symbols in terms of components of the metric 
tensor are 

»B 2 9 V d& d£ a d£* 
The following equality takes place: 



a 1 dJg 

r^ = r^ = — — , ^ = det||^ 



The covariant derivative of a scalar function coincides with partial derivatives 
with respect to the corresponding coordinate £° : . The covariant derivative of 
the contravariant coordinate of a vector field has the form 



Appendix B. Group classification 



Here the details of the group classification of system ([15]) with respect 
to function h(p, p) as an "arbitrary element" are given. Results of the ad- 
missible group calculations are summed up in Theorem [2] The classifying 



equation is (17). One need to specify all possible forms of function h(p, p) 



up to equivalence transformations (16), and the corresponding restrictions 
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on arbitrary constants q such that equation (17) is satisfied identically for p 
and p. 

The classifying parameter is the dimension dim{V} of the linear space 
spanned by the vector 

V = (h + Aph p + 5ph p , ph p + ph p , h p ) 

for various values of p and p. 
a) dim{V} = 0. This implies 



h + Aph p + 5ph p = 0, ph p + ph p = 0, h f 



0. 



Hence, function h is identically zero. Constants c%, C2, C3 take arbitrary 
values. This case gives the maximal extension of the admissible group, 
b) dim{V} = 1. Hence, vector V is proportional to the constant vector: 

h + 4ph p + 5ph p = mU(p, p), ph p + ph p = nU(p, p), h p = sU(p,p). 

The equivalent form of these equations is 

h = (m — 5n + ps)U, h p = p~ l {n — ps)U, h p = sU. 

Calculation of mixed derivatives yields the following equations for function 

U(p,p): 



(m — 5n + ps)U p = 0, p 1 (n — ps)U = (m — 5n + ps)U l 
sU + spU p — (n — ps)U p . 



(B.l) 



First of equations (B.l) and condition h ^ give U p = 0, hence 

U = U{p). 

To satisfy the last equation of ( B.l[ ) the following two cases must be sepa- 
rated: s = and s ^ 0. 

s = . Solution of the remaining second equation in (B.l) is 

n 

U = p m - 5n Uo. 



Accurate to the equivalence transformations (16) function h can be brought 
to the form 

h = p k . 
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Hence, vector V is equal to 



V = p k (l + 5fc,fc,0). 



The classifying equation (17) is reduced to 



Ci(l + 5k) — 2kc 2 = 0, c 3 is arbitrary. 
Hence, the admissible group is extended by two operators. 



s 7^ 0. The last of equations (B.l ) gives 

U 



U 



P 



From the second of equations (B.l) it follows that m = An. Accurate to the 



equivalence, function h can be brought to the form 

h = kp/p. 
with an arbitrary constant k. Hence, 

V = (0,0,l/p) 

which implies that C\ and c 2 are arbitrary constants, and C3 = 0. The group 
is extended by two operators. 

c) dim{V} = 2. In this case there exists a constant vector (to, n, s) 
orthogonal to vector V for every p and p, that is 

m(h + Aph p + 5ph p ) + n(ph p + ph p ) + sh p = 0. (B.2) 

Integration of this equation in the case 4m + n^0 for function h(p, p) gives 

— m / 5m-\-n 

h = ((4m + n)p + s) 4m+n / f p ((4m + n)p + s) 4m+n 



with an arbitrary function /. Accurate to transformations (16) this is equiv- 
alent to 

h = p k f(pp k ~ 1 ). (B.3) 



with some constant parameter k. Equation (17) is satisfied only for 



1 + 4k)a - 2c 2 k = 0, c 3 = 0. 
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The admissible group is extended by one operator. In the case Am + n = 0, 
s the general solution of equation (B.2) has the following form: 

h = e- mp/s f(pe- mp/s ), 

which is equivalent to either 

h = f(pe p )/ P 

for m 7^ 0, or 



h = tip) 



(B.4) 
(B.5) 



otherwise. Equation (17) gives 



for the case (B.4), and 



C2 — C3 — 2ci 



c\ = c 2 = 



for the case (B.5). Finally, the case 4m + n = 0, s = gives the following 
solution 

h = f(p)/p. (B.6) 



The classifying equation (17) yields 

c 2 = 2ci, c 3 = 0. 

Hence, class of function h which gives two-dimensional vector space {V} up 
to the equivalence is exhausted by representatives (B.3)-(B.6). 

d) dim{V} = 3. This case is equivalent to an arbitrary function h. The 
admissible group coincides with the kernel. 



References 

[1] H. K. Moffatt. Field Generation in Electrically Conducting Fluids. Cam- 
bridge University Press, Cambridge, London, New York, Melbourne, 
1978. 

[2] V. I. Arnold and B. A. Khesin. Topological methods in hydrodynamics, 
v. 125 of Applied Mathematical Sciences. Springer- Verlag: New York, 
1998. 



23 



[3] M. A. Berger, L. H. Kauffman, B. Khesin, H. K. Moffatt, and R. L. 
Ricca. Lectures on topological fluid mechanics: lectures given at the 
CIME Summer School held in Cetraro, Italy, July 2-10, 2001. Springer 
Verlag, 2009. 

[4] S. V. Golovin. Analytical description of stationary ideal MHD flows 
with constant total pressure. Phys. Lett. A., 374:901-905, 2010. 

[5] C. Rogers and W. K. Schief. Novel integrable reductions in nonlin- 
ear continuum mechanics via geometric constraints. J. Math. Phys., 
44(8):3341-3369, 2003. 

[6] D. Gilbarg. On the flow patterns common to certain classes of plane 
fluid motions. J. Math. Phys., 26:137, 1947. 

[7] W. K. Schief. Hidden integrability in ideal magnetohydrodynamics: The 
Pohlmeyer-Lund-Regge model. Physics of plasmas, 10(7):2677-2685, 
2003. 

[8] J. D. Gibbon and D. D. Holm. Lagrangian analysis of alignment dynam- 
ics for isentropic compressible magnetohydrodynamics. New Journal of 
Physics, 9:292, 2007. 

[9] B. F. Schutz. Geometrical methods of mathematical physics. Cambridge: 
Cambridge University Press, 1980. 

[10] L. V. Ovsiannikov. Group analysis of differential equations. Academic 
Press: New York, 1982. 

[11] P. J. Olver. Applications of Lie groups to differential equations. Springer 
Verlag, 2000. 

[12] S. Chandrasekhar. On the stability of the simplest solution of the equa- 
tions of hydromagnetics. Proc. Nat. Acad. Set. U.S.A., 42:273-276, 1956. 



24 



